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摘 要 : 掌握 陆地 水 储量 变化 的 空间 分 布 特点 和 变化 情况 ,了 解 其 变化 的 原因 有 助 于 水 资源 的 可 持 
续 性 综合 管理 。 以 中 国 陆 地 为 研究 区 域 , 基 于 重力 恢复 和 气候 实验 (GRACE) 卫 星 反 演 得 到 的 陆地 
水 储量 变化 数据 ,使 用 曼 - 肯 德尔 法 (Mann-Kendall, M-K) 趋 势 检 验 和 经 验 正 交 函数 (EOF) 分 析 陆 


地 水 储量 变化 的 趋势 以 及 时 空 变化 特征 ,并 选择 了 包括 气象 下 执 


和 社会 经 济 3 个 方面 的 10 个 


影响 因子 , 即 气 温 、 降 水 标准 化 降水 蒸 散发 指数 (SPEI) 、 不 透水 层 占 比 、 水 体 占 比 、 归 一 化 植被 指 
ZCONDVD HE 坡度 、 国 内 生产 总 值 和 人 口 ,利用 地 理 探 测 器 、 皮 尔 逊 (Pearson) 相 关 分 析 和 随机 
森林 3 种 方法 综合 分 析 其 与 陆地 水 储量 变化 的 关系 。 结 果 表 明 : 松 花 江 、 嫩 江 和 松 嫩 平 原 附 近 以 及 


柴 达 木 盆 地 -长 江 -东南 没 海 条 带 上 陆地 水 储量 


( 即 0 个 月 滞后 ) 为 主 。 


变化 显著 增加 ,中 国 西南 以 及 新 疆 -黄土 高 原 -华北 
平原 条 带 上 陆地 水 储量 变化 显著 减少 ,从 高 续 度 到 低 纬度 呈现 出 “高 - 低 -高 - 低 " 的 变化 特点 。 整 
体 而 言 , 气 象 因素 对 陆地 水 储量 变化 的 解释 力 最 强 ,其 次 是 社会 经 济 因素 和 下 垫 
量 是 引起 陆地 水 储量 变化 的 主要 因素 。 月 陆地 水 储量 变化 对 降水 量 ` 温 度 SPEI 和 NDVI 响 应 存在 
滞后 性 ,上 且 陆 地 水 储量 变化 对 降水 量 和 SPEI 以 1 个 月 沾 后 为 主 ,对 温度 和 NDVI 的 响应 以 无 滞后 


因素 ,其 中 降水 


X 键 词 : 陆地 水 储量 变化 ; CRACE 卫星 ; 影响 分 析 ; 滞后 效应 
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陆地 水 储量 通常 包括 地 下 水 .土壤 水 .地表 水 
(湖泊 河流 和 水 库 ) 冰川 与 积 雪 等 部 分 上 。 陆 地 水 
储量 的 变化 通常 被 认为 是 一 个 关键 的 水 文 指标 , 它 
反映 了 整个 地 区 水 文 输入 (如 降雨 降雪、 冰川 融化 
和 径流 ) 与 输出 (如 荧 散 、 人 类 用 水 和 径流 ) 之 间 的 
净 效应 ,影响 下 渗 与 基 流 过 程 。 掌 握 陆 地 水 储量 变 
化 的 时 空 特 征 , 了 解 其 变化 的 原因 有 助 于 水 资源 的 
可 持续 性 综合 管理 。 

陆地 水 储量 变化 的 增 减 会 改变 水 文 循环 ,威胁 
生态 系统 与 人 口 的 可 持续 用 水 。 美 国 国家 航空 航 
天 局 (NASA ) 和 德国 航空 航天 中 心 (DLR) 于 2002 年 
3 月 联合 开发 了 重力 恢复 和 和 气候 实验 (GRACE ) 卫 
星 , 利 用 GRACE 卫星 观察 到 的 地 球 重 力 场 变化 ,可 
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反 演 地 球 系统 中 的 陆地 水 储量 变化 ,这 大 大 方便 了 
对 陆地 水 储量 变化 的 研究 中 。 迄 今 为 止 ,不 少 学 者 
针对 陆地 水 储量 变化 的 时 空 特征 和 影响 因素 开展 
了 研究。 相关 研究 指出 中 国 陆地 地 区 北方 陆地 水 
储量 变化 以 减少 为 主 ,南方 以 增加 为 主 ”” ,云南 省 、 
天 山 山区 以 及 西北 地 区 陆地 水 储量 变化 整体 上 呈 
现下 降 趋势 , 目 天 山 山区 以 及 西北 地 区 的 陆地 水 储 
量变 化 呈现 的 明显 季节 变化 主要 由 降水 造成 "1。 
Yang 等 中 Wei 等 "分 别 对 塔里木 河和 柴 达 木 伟 地 陆 
地 水 储量 进行 了 影响 因素 分 析 , 均 指出 降水 变化 是 
所 研究 区 域 陆 地 水 储量 变化 的 主要 驱动 因素 。 李 
晓 英 等 中 从 趋势 性 、 相 关 性 及 时 空 变化 特征 3 个 方 
面 分 析 了 长 江 流 域 陆地 水 储量 变化 以 及 与 植被 履 


作者 简介 : 石 振 君 (1998-) , 女 ,硕士 研究 生 ,主要 从 事 遥 感应 用 研究 . E-mail: 16678134329@163.com 
通讯 作者 : 朱 秀 芳 (1982- ) , 女 ,博士 ,教授 ,主要 从 事 遥 感应 用 研究 . E-mail: zhuxiufangGbnu.edu.cn 


202309.00233v1 


chinaXiv 


ChinaXiv 合 作 期 刊 


1398 T 4 w 4635 


盖 变 化 的 关系 ,指出 长 江 流域 归 一 化 植被 指数 与 陆 
地 水 储量 变化 之 间 有 强 相关 性 。Meng 等 "从 地 表 水 
平衡 的 角度 出 发 ,研究 发 现 青 藏 高 原 陆地 水 储量 的 
变化 主要 归 因 于 降水 和 莹 散 的 变化 。Zhong 等 王 在 
海河 流域 的 研究 表明 每 年 由 气候 驱动 的 陆地 水 储 
量 虽 呈现 一 定 的 波动 上 升 趋势 ,但 长 期 的 人 为 取水 
导致 陆地 水 储量 持续 下 降 。 

以 上 研究 大 多 在 中 国 的 局 部 区 域 开展 ,还 没有 
看 到 对 整个 中 国 地 区 的 陆地 水 储量 变化 的 影响 因 
素 分 析 。 为 此 ,本 文 使 用 GRACE 卫星 反 演 得 到 的 
陆地 水 储量 变化 数据 ,以 中 国 陆地 为 研究 区 域 , 探 
究 近 十 几 年 来 陆地 水 储量 变化 ,综合 分 析 气 象 因 
素 、 下 热 面 和 社会 因素 对 陆地 水 储量 变化 的 影响 ， 
为 我 国 水 资源 管理 提供 适用 的 评估 工具 和 科学 决 
策 支 撑 信 息 。 


1 数据 与 方法 


1.1 研究 区 概况 

中 国 位 于 亚洲 的 东部 ,太平 详 的 西岸 ,地 势 呈 
西高 东 低 的 阶梯 状 分 布 。 中 国 经 纬度 跨度 大 ,受到 
气候 .地势 等 因素 的 共同 影响 ,中 国 的 水 资源 分 布 
呈现 出 东 多 西 少 . 南 多 北 少 的 格局 。 中 国 的 东部 和 
南部 地 区 主要 是 亚热带 季风 和 温带 季风 气候 ,降水 
呈现 出 明显 的 季节 性 特征 ,夏季 和 秋季 降水 多 , 春 
季 和 冬季 降水 少 ;西部 和 北部 地 区 主要 为 温带 大 陆 
性 气候 和 高 原 山 地 气候 ,全 年 降水 均 比 较 少 。 
1.2 数据 来 源 

本 文 所 用 到 的 数据 包括 陆地 水 储量 变化 数据 、 
气象 数据 地表 数 据 和 社会 经 济 数据 。 

陆地 水 储量 变化 数据 "来 自 国家 青藏 高 原 科 
学 数据 中 心 ,为 2003 一 2019 年 的 中 国 区 域 陆地 水 储 
量变 化 月 数据 ,空间 分 辩 率 为 0.23"。 该 数据 集 连 接 
了 CSR RL06 Mascon 和 JPL RL06 Mascon 数据 在 
GRACE 和 GRACE-FO 之 间 的 11 个 月 的 间断 期 ,但 
并 未 对 GRACE 和 GRACE-FO 各 自 数据 之 间 存 在 的 
个 别 月 份 的 缺失 数据 进行 弥补 。 在 使 用 数据 前 ,使 
用 线性 插值 方法 "对 个 别 月 份 的 缺失 数据 进行 了 
插 补 ,得 到 了 完整 的 陆地 水 储量 变化 时 序数 据 。 

气象 数据 包括 降水 量 温度 和 标准 化 降水 芋 散 
发 指数 (SPET) 。 降 水 量 和 温度 数据 均 为 来 自 国家 
气象 数据 中 心 的 2003 一 2019 年 的 月 数据 ,空间 分 辩 


率 为 0.23"。SPEI 数 据 为 来 自 全 球 标准 化 降水 蒸 散 
指数 数据 库 的 2003 一 2019 年 的 月 数据 ,空间 分 辩 率 
为 0.5°。 

地 表 数 据 包括 土地 利用 现状 图 、 归 一 化 植被 指 
数 (NDVI) 和 数字 高 程 模型 (DEM) 数 据 。 土 地 利用 
现状 图 为 来 自 于 资源 环境 科学 与 数据 中 心 的 1 km 
分 辩 率 的 数据 ,时 间 分 别 为 2005 .2010 年 和 2015 
年 。NDVI 数 据 为 NASA 的 陆地 过 程 分 布 式 数据 档 
案 中 心 提供 的 2003 一 2019 年 的 MODIS 的 陆地 标准 
产品 ,空间 分 辩 率 为 1 km ,时 间 分 辨 率 为 8 d。NDVI 
数据 包括 年 度 NDVI 数 据 集 和 月 度 NDVI 数 据 集 。 
月 度 NDVI 数 据 集 通过 最 大 值 合成 法 生成 ,年 度 ND- 
VI 数据 集 是 在 月 度数 据 基础 上 采用 平均 值 合成 法 
生成 。DEM 为 基于 最 新 的 SRTM V4.1 数据 经 重 采 
样 生成 的 分 辩 率 为 1 km 的 数据 。 社 会 经 济 数据 包 
括 国 内 生产 总 值 (GDP) 数 据 和 人 口 数 据 ,2 类 数据 
均 来 自 与 资源 环境 科学 与 数据 中 心 ,分 辩 率 均 为 
1 km, 所 用 数据 的 时 间 分 别 为 2005、2010 年 和 2015 
年 。 在 使 用 前 ,对 上 述 数 据 进行 了 重 投影 和 重 采样 
处 理 , 重 采样 后 的 分 辨 率 为 0.25°。 
1.3 研究 方法 

(1) 曼 - 肯 德尔 法 (Mann-Kendall, M-K ) 趋 热 检验 

M-K 趋势 检验 是 一 种 非 参 数 方法 ,其 目的 是 
评估 变量 随 着 时 间 的 变化 是 否 有 单调 上 升 或 下 降 
H5 30 。M-K 检 验 不 需要 样本 遵循 特定 的 分 
布 ,也 不 会 受 少数 异常 值 的 干扰 ,即使 数据 有 缺失 
值 或 者 存在 值 低 于 一 个 或 者 多 个 检测 的 限制 ,也 可 
以 使 用 M-K 检 验 。 

(2) 经 验 正 交 也 数 (EOF) 分 析 方 法 

EOF 分 析 方 法 是 提取 数据 主要 特征 量 的 一 种 
PU, EOF 分 解 得 到 的 特征 向 量 对 应 的 是 空间 
样本 ,也 称 作 空 间 特征 向 量 或 者 空间 模 态 ,反映 地 
理 要 素 的 空间 分 布 特点 ; 主 成 分 对 应 的 是 时 间 变 
化 ,也 称 时 间 系 数 , 反 映 相 应 空间 模 态 随时 间 的 变 
化 而 发 生 的 变化 。 

(3) 皮尔 逊 (Pearson) 相 关系 数 

Pearson 相关 系数 是 一 种 线性 相关 系数 ,也 是 最 
篆 用 的 一 种 相关 系数 拉 。 其 值 介 于 -1 到 1 之 间 , 绝 
对 值 越 大 表明 相关 性 越 强 。 值 大 于 0 时 ,表明 2 个 
变量 正 相关 ; 值 小 于 0 时 ,表明 2 个 变量 负 相 关 ; 当 
等 于 0 时 ,表明 2 个 变量 不 是 线性 相关 的 ,但 是 可 能 
存在 其 他 方式 的 相关 。 


202309.00233v1 


chinaXiv 


(4) 随机 森林 

随机 森林 是 基于 Bagging 框 架 通 过 集成 学 习 的 
思想 将 多 棵 树 进 行 集成 的 一 种 算法 ,其 基本 单元 是 
决策 树 ,本 质 是 集成 学 习 方 法 。 随 机 森林 模型 不 易 
发 生 过 拟 合 问题 ,具有 较 强 的 泛 化 能 力 。 它 可 以 用 
于 评估 各 个 特征 对 于 结果 的 重要 性 2 。 

(5) 地 理 探测 分 析 

地 理 探测 器 是 探测 因子 的 空间 分 异性 ,并 揭示 
其 驱动 力 的 一 组 方法 中 。 地 理 探测 器 的 核心 思想 
是 如 果 某 个 自 变 量 对 某 个 因 变量 有 重要 的 影响 , 那 
么 该 自 变 量 和 因 变 量 的 空间 分 布 应 该 具有 相似 
性 。 地 理 探测 器 既 可 以 用 于 数值 型 数据 ,也 可 以 用 
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于 定性 数据 。 地 理 探 测 器 通过 计算 和 比较 各 单 因 
子 的 g 值 以 及 两 因子 受 加 后 的 g 值 ,判断 两 因子 是 否 
存在 交互 作用 ,以 及 交互 作用 的 强 弱 方向、 线性 还 
是 非 线 性 等 。 地 理 探测 器 包括 分 异 及 因子 探测 、 交 
互 作用 探测 .风险 区 探测 和 生态 探测 4 个 子 探测 需 。 
1.4 技术 路 线 

本 文 的 技术 路 线 (图 1) 主 要 包括 3 个 内 容 : 基 于 
M-K 趋势 检验 的 陆地 水 储量 变化 的 趋势 分 析 ,基于 
EOF 的 陆地 水 储量 变化 的 时 空 特征 分 析 , 基 于 Pear- 
son 相关 、 随 机 森林 和 地 理 探测 器 的 陆地 水 储量 变 
化 影响 因素 分 析 。 其 中 ,陆地 水 储量 变化 的 影响 
素 分 析 的 主要 步骤 如 下 : 


影响 因子 
土地 利用 现状 图 国内 生产 总 
气象 数据 MODIS NDVI 
SRTM DEM AH 
GRACE 陆地 水 储量 y Y 
变化 数据 
下 执 面 因素 
气象 因素 (不 透水 层 占 比 、 水 社会 经 济 因素 
(气温 、 降 水 量 和 SPED 体 占 比 、 m 高 程 (GDP 和 人 口 ) 
和 坡度 ) 


Mann-Kendall Ze EAS PRÉC 、 


变化 趋势 时 空 变化 特征 影响 因素 重要 性 排序 和 时 间 滞 后 效应 


注 :SPEI 为 标准 化 降水 莹 散发 指数 ;NDVI 为 归 一 化 植被 指数 ;DEM 为 数字 高 程 模型 ;CDP 为 国内 生产 总 值 。 


图 1 技术 路 线 


Fig. 1 Technical route 


(1) 影响 因子 选择 :包括 气象 因素 (气温 、 降 水、 
SPEL) \ 下 热 面 因 素 ( 不 透水 层 占 比 水体 占 比 、ND- 
VI ,高程 .坡度 ) 和 社会 经 济 因素 (CDP、 人 口 )3 个 方 
面 10 个 影响 因子 。 

(2) Pearson 相关 分 析 : 在 10 个 影响 因子 中 ,ND- 
VI, SPEI, 降水 量 和 温度 会 随 着 时 间 的 变化 而 发 生 
较 大 变化 , 且 存 在 时 间 序 列 的 数据 ,其 他 6 个 因子 在 
短 时 间 内 的 变化 量 很 小 , 且 没 有 时 间 序 列 的 数据 。 
因此 ,基于 了 Pearson 的 相关 分 析 分 成 了 2 部 分 ,一 部 
分 是 时 间 维 的 分 析 , 另 一 部 分 是 空间 维 的 分 析 。 时 
间 维 的 分 析 是 考虑 NDVI、SPEI、 降 水 量 .温度 与 陆 
地 水 储量 变化 存在 滞后 效应 , 逐 栅 格 计算 NDVI、 
SPEI, 降水 量 温度 与 陆地 水 储量 变化 在 0~3 个 月 滞 


后 时 间 内 的 相关 系数 ,并 找 出 最 大 清 相 关系 数 对 应 
的 滞后 月 数 。 空 间 维 的 分 析 是 首先 过 机 格 计算 时 
间 序 列 数据 (包括 NDVI.SPEI、 降 水 量 .温度 和 陆地 
水 储量 变化 ) 的 多 年 均值 ,然后 以 所 有 格 网 为 样本 ， 
分 别 计算 陆地 水 储量 变化 和 10 个 影响 因子 的 Pear- 
son 相关 系数 ,以 相关 系数 大 小 来 反映 影响 因子 和 
陆地 水 储量 变化 的 相关 程度 。 

(3) 基于 随机 森林 进行 影响 因子 重要 性 排序 : 
对 随机 森林 中 某 个 树 节 点 的 样本 方差 与 分 裂 后 2 个 
叶子 节点 样本 方差 进行 差 值 运算 ,用 该 差 值 来 评价 
某 个 特征 的 重要 性 程度 , 差 值 越 高 表示 该 特征 重要 
性 越 大 。 计 算得 到 模型 所 有 特征 重要 性 后 ,对 所 有 
特征 作 归 一 化 处 理 , 对 于 某 一 模型 ,所 有 特征 重要 
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性 值 累计 为 1 ,特征 重要 性 具体 的 值 本 身 并 无 意义 ， 
因此 分 析 时 关注 的 是 特征 重要 性 之 间 的 相对 差异 。 

(4) 基于 地 理 探测 噩 探 究 单个 影响 因子 及 因子 
间 的 交互 作用 对 陆地 水 储量 变化 空间 差异 的 解释 
71 :在 使 用 地 理 探测 絮 计 算 之 前 ,和 完 对 自 变量 进行 
离散 化 处 理 ,按照 干旱 等 级 划分 标准 对 平均 SPEI 进 
行 分 级 ,对 其 他 9 个 影响 因子 (年 均 NDVI\ 年 均 降 水 
量 \ 年 均 温度 、 人 口 .GDP 不 透水 层 占 比 水 体 占 比 ， 
高 程 和 坡度 ) 用 自然 间断 点 分 级 法 重 分 类 为 6 类 。 
然后 以 离散 化 后 的 因子 为 *, 以 地 表 水 储量 变化 趋 
势 值 为 ,进行 因子 探测 .交互 探测 和 生态 探测 分 
析 。 因 子 探测 用 来 确定 10 个 影响 因子 多 大 程度 上 
解释 了 陆地 水 储量 变化 的 空间 分 异 ; 交 互 探测 用 于 
评估 2 个 自 变量 共同 作用 时 是 否 会 增加 或 减弱 对 因 
变量 的 解释 力 ;生态 探测 用 于 比较 2 个 影响 因子 之 
间 对 陆地 水 储量 变化 的 空间 分 布 的 影响 是 否 有 显 
车 的 差异 。 


2 结果 与 分 析 


2.1 陆地 水 储量 变化 趋势 

由 2003 一 2019 年 中 国 陆地 水 储量 变化 趋势 图 
(图 2) 可 以 看 出 ,大 部 分 地 区 陆地 水 储量 变化 趋势 
明显 ,其 中 陆地 水 储量 变化 显著 增加 的 区 域 主要 分 
布 在 松花 江 嫩江 和 松 嫩 平原 附近 以 及 柴 达 木 盆 
地 -长 江 - 东 南 沿海 条 带 上 ,陆地 水 储量 变化 显著 
减少 的 区 域 主要 分 布 在 中 国 西南 以 及 新 疆 - 黄 土 


m 显著 减少 
国 不 显著 
tA VUE m 显著 增加 
审 图 号 : GS(2023)2267 号 


0 600 km 
L_ 


图 2 中 国 陆地 水 储量 变化 趋势 


Fig.2 Change trend of terrestrial water storage in China 


高 原 - 华 北平 原 条 带 上 。 
2.2 陆地 水 储量 变化 时 空 特征 

对 2003 一 2019 年 的 204 幅 (每 月 1 幅 影 像 ) 陆 地 
水 储量 变化 影像 进行 EOF 分 析 。EOF 的 前 3 个 模 态 
的 解释 方差 分 别 为 47.98%、19.14% ,16.8996 , RIFI 
差 贡 献 约 为 84.01%。 由 图 3a~c 可 以 看 出 ,第 1 模 态 
表现 出 中 国 的 陆地 水 储量 变化 随 纬度 带 交 替 反 相 
变化 的 特征 ,从 高 纬 到 低 纬 呈 现 出 水 储量 变化 “高 - 
低 - 高 - 低 ” 的 变化 特征 ;第 2 模 态 表现 出 中 国 南 北方 
陆地 水 储量 变化 反 相 变化 的 特征 ;第 3 模 态 显示 出 
四 川 盆 地 横断 山脉 .云贵 高 原 、. 某 达 木 盆地 陆地 水 
储量 变化 比较 高 ,而 新 疆 部 分 区 域 和 中 国 东南 部 水 
储量 变化 明显 偏 低 。 结 合 时 间 系 数 看 (图 3d~f) ,第 
1 模 态 的 时 间 系 数 在 研究 时 间 段 的 前 9 a 均 为 负 值 ， 
而 后 6 a 均 为 正 值 ,数值 呈现 逐渐 增 大 的 特征 ,说明 
随 着 时 间 的 变化 ,中 国 陆地 的 水 储量 变化 由 第 1 模 
态 反 应 的 由 北 到 南 的 “高 - 低 - 高 - 低 ” 特 征 的 反 向 逐 
渐 演 变 到 由 北 到 南 的 “高 - 低 - 高 - 低 ” 的 空间 格局 特 
征 。 第 2 模 态 和 第 3 模 态 反应 出 的 特征 具有 震荡 变 
化 的 特点 ,总体 来 说 每 年 的 上 半年 的 时 间 系 数 为 负 
值 , 下 半年 的 时 间 系 数 为 正 值 ,说 明 上 半年 的 各 月 
的 陆地 水 储量 与 第 3 模 态 反应 的 特征 相反 ,而 下 半 
年 与 第 3 模 态 反应 的 特征 一 致 。 
2.3 影响 因素 分 析 
2.3.1 基于 地 理 探 测 器 的 影响 因子 重要 性 分 析 
子 探测 .交互 探测 和 生态 探测 的 结果 如 表 1~3 所 
示 。 在 单个 因子 中 ,降水 量 对 陆地 水 储量 变化 解释 
力 最 强 ,g 值 为 0.276; 水 体 占 比 对 陆地 水 储量 变化 
解释 力 最 弱 ,g 值 仅 为 0.001。 在 因子 交互 作用 中 ， 
对 陆地 水 储量 变化 解释 力 最 强 的 是 降水 量 和 温度 
的 交互 作用 ,g 值 达到 0.461; 其 次 是 降水 量 和 DEM 
的 交互 作用 ,g 值 为 0.377; 水 体 占 比 和 坡度 的 交互 
作用 对 陆地 水 储量 变化 的 解释 力 最 弱 ,g 值 为 
0.037。 综 合 比 较 因子 探测 和 交互 探测 的 结果 ,可 
得 到 因子 两 两 交互 作用 相 较 于 单个 因子 对 于 陆地 
水 储量 变化 的 解释 力 均 为 增加 。 生 态 探 测 结 果 表 
HH ,降水 量 与 DEM .GDP、NDVI、 人 口 ,SPEI 与 GDP、 
AH , ii EE 5j DEM .GDP、NDVI、 人 口 .SPEI 对 陆地 
水 储量 变化 空间 分 布 的 影响 有 显著 的 差异 ,其 余 2 
个 因子 之 间 对 陆地 水 储量 变化 空间 分 布 的 影响 无 
显著 差异 。 
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图 3 中 国 陆地 水 储量 变化 的 EOF 分 析 
Fig.3 EOF analysis of terrestrial water storage change in China 
表 1 因子 探测 结果 
Tab.1 Factor detection results 
影响 因子 q HEF 影响 因子 4 值 排序 
降水 量 0.276 1 DEM 0.042 6 
温度 0.149 2 AH 0.035 7 
SPEI 0.064 3 GDP 0.028 8 
NDVI 0.050 4 坡度 0.017 9 
不 透水 层 占 比 0.044 5 水 体 占 比 0.001 10 


注 :SPEI 为 标准 化 降水 蒸 散 发 指数 ;NDVI 为 归 一 化 植被 指数 ;DEM 为 数字 高 程 模型 ;GDP 为 国内 生产 总 值 ;4 值 为 影响 因子 解释 力 。 下 同 。 
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表 2 因子 交互 探测 结果 
Tab.2 Factor interactive detection results 

影响 因子 DEM GDP NDVI AH 降水 量 SPEI 温度 不 透水 层 占 比 KRE ”坡度 
DEM 0.042 - - - - - - - - - 
GDP 0.085 0.028 - - - - - - 一 一 
NDVI 0.116 0.092 0.050 - - - - - - - 
A 0.096 0.049 0.102 0.035 - - - - - - 
降水 量 0.377 0.343 0.349 0.353 0.276 = - - - - 
SPEI 0.184 0.170 0.146 0.172 0.373 0.064 - - - 一 
温度 0.308 0.218 0.220 0.233 0.461 0.341 0.149 - - - 
不 透水 层 占 比 0.140 0.082 0.117 0.093 0.337 0.143 0.248 0.044 一 一 
水 体 占 比 0.069 0.043 0.059 0.053 0.308 0.082 0.186 0.075 0.001 一 
坡度 0.093 0.056 0.070 0.063 0.293 0.087 0.190 0.068 0.037 0.017 

表 3 生态 探测 结果 
Tab. 3 Ecological detection results 

影响 因子 DEM GDP NDVI 人 降水 量 SPEI 温度 不 透水 层 占 比 。 坡度 
GDP - - - - - - - 
NDVI N - - - - - - - 
AH N - - - - - 
降水 量 Y Y Y Y - - - - - 
SPEI Y Y N - - - - 
温度 Y Y Y Y N Y - - - 
不 透水 层 占 比 N N N N N - = 
坡度 N N N N N N N N - 
水 体 占 比 N N N N N N N 


TE Y N 分 别 表示 2 个 影响 因子 对 陆地 水 储量 变化 空间 分 布 的 影响 丰 


2.3.2 基于 随机 森林 的 影响 因子 重要 性 分 析 由 基 


E95% 显 著 性 水 平 上 存在 和 不 存在 明显 差异 。 


于 陆地 水 储量 变化 的 影响 最 大 ,其 次 是 社会 经 济 


于 随机 和 森林 的 影响 因子 重要 性 分 析 结 果 ( 表 4) 可 
知 ,10 个 影响 因子 对 陆地 水 储量 变化 的 重要 性 排序 
为 降水 量 > 温 度 >SPEI> 人 口 >GCDP>DEM>NDVI> 坡 
度 > 不 透水 层 占 比 > 水 体 占 比 ,降水 量 因子 对 陆地 水 
储量 变化 的 影响 最 大 ,特征 重要 性 值 为 0.335 ,水 体 
占 比 因子 对 陆地 水 储量 变化 的 影响 最 小 ,特征 重要 
性 值 仅 为 0.032。 在 因子 类 型 层面 分 析 , 气 象 因 子 对 


子 , 下 热 面 因子 对 于 陆地 水 储量 变化 的 影响 最 小 。 

2.3.3 Æ F Pearson 相关 分 析 的 影响 因子 重要 人 性 分 
析 由 空间 维 上 陆地 水 储量 变化 和 各 影响 因子 相 
关 性 分 析 的 结果 ( 表 5) 可 知 ,从 全 区 来 看 ,DEM、 人 
口 和 水 体 占 比 因子 与 陆地 水 储量 变化 不 存在 线性 
相关 关系 ;不 透水 层 .坡度 .GDP 与 陆地 水 储量 变化 
呈 负 相关 , 且 不 透水 层 与 陆地 水 储量 变化 的 负 相关 


表 4 基于 随机 森林 的 影响 因子 重要 性 评价 结果 


Tab.4 Evaluation results of importance of influencing factors based on random forests 
影响 因子 重要 性 值 排序 影响 因子 重要 性 值 排序 
降水 量 0.335 1 DEM 0.072 6 
温度 0.160 2 NDVI 0.059 7 
SPEI 0.110 3 坡度 0.039 8 
人 0.082 4 不 透水 层 占 比 0.038 9 
GDP 0.073 5 水 体 占 比 0.032 10 
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X5 Pearson 相关 分 析 结 


Tab. 5 Pearson correlation analysis results 


影响 因子 一 一 GOES 影响 因子 
相关 系数 显著 性 检验 排序 相关 系数 显著 性 检验 排序 
降水 量 0.339 是 1 GDP -0.035 是 6 
SPEI 0.207 是 2 温度 0.030 是 7 
不 透水 层 -0.160 是 3 DEM -0.014 fr 8 
坡度 -0.101 是 4 AH -0.012 fr 9 
NDVI 0.048 是 5 水 体 占 比 -0.005 fi 10 


注 : 表 中 排序 是 按照 相关 系数 的 绝对 值 进行 的 排序 。 


性 最 强 ,相关 系数 为 -0.16;NDVI\ 降 水 量 .SPEI 和 温 
度 与 陆地 水 储量 变化 呈正 相关 关系 , 且 降 水 量 与 陆 
地 水 储量 变化 的 正 相 关 性 最 强 ,相关 系数 为 0.339。 

陆地 水 储量 变化 对 气候 因素 和 NDVI 的 响应 沛 
后 时 间 分 布 如 图 4 所 示 。 大 部 分 区 域 陆 地 水 储量 变 
化 与 NDVI 在 当月 以 及 滞后 1 个 月 时 相关 性 最 强 , 陆 
地 水 储量 变化 与 NDVI 存 在 2~3 个 月 沸 后 的 区 域 主 
要 位 于 中 国 中 部 地 区 (图 4a)。 中 国 大 部 分 地 区 的 
陆地 水 储量 变化 与 降水 量 存在 1 个 月 的 滞后 (图 


(a) 陆地 水 储量 变化 与 NDVI 滞 相关 人 


图 例 一 一 国界 v 未 定 国 界 
滞后 时 间 EH 无 沸 后 


4b)。 具 体 来 说 ,中 国 西部 靠 北 地 区 陆地 水 储量 变 
化 与 当月 降水 量 相关 性 最 高 ,中 国 东 北部 、 东 南部 
以 及 西部 部 分 地 区 陆地 水 储量 变化 在 滞后 1 个 月 时 
与 降水 量 相关 性 最 高 ,中 国 中 部 地 区 陆地 水 储量 变 
化 在 滞后 2 个 月 时 与 降水 量 相关 性 最 高 ,华北 平原 
中 部 河北 以 及 河北 北部 地 区 陆地 水 储量 变化 在 在 
滞后 2 个 月 时 与 降水 量 相关 性 最 高 。 中 国 大 部 分 地 
区 的 陆地 水 储量 变化 与 温度 存在 0 或 3 个 月 的 沾 后 
(图 4c)。 具 体 来 说 ,中 国 东北 部 、 东 南部 以 及 西北 


(6) 陆地 水 储量 变化 与 降水 量 灌 相 关 À 


省 界 台湾 省 资料 暂 缺 
滞后 1 个 月 国 一 滞后 2 个 月 HB 滞后 3 个 月 


审 图 号 : GS(2023)2267 号 


图 4 陆地 水 储量 变化 与 影响 因子 之 间 的 滞 相 关 分 析 结 果 


Fig.4 Lag correlation analysis results between terrestrial water storage change and influencing factors 
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部 地 区 陆地 水 储量 变化 与 当月 的 温度 相关 性 最 高 ， 
青藏 高 原 中 部 和 南部 陆地 水 储量 变化 在 滞后 2 个 月 
时 与 温度 相关 性 最 高 , 云南 省 、 华 北平 原 以 及 华北 
平原 西部 地 区 陆地 水 储量 变化 在 沾 后 3 个 月 时 与 温 


降水 量 是 导致 陆地 水 储量 变化 的 主要 原因 ,这 与 
Wei 等 ”Meng 等 "研究 结论 一 致 。 但 需 注意 的 是 ， 
降水 量 增加 并 不 能 完全 代表 陆地 水 储量 增加 ,例如 
新 疆 部 分 地 区 降水 量 增 加 明显 ,但 蒸 散发 能 力 明 显 


度 相关 性 最 高 。 中 国 大 部 分 地 区 的 陆地 水 储量 变 
化 与 SPEI 存 在 1 或 3 个 月 的 滞后 (图 4d)。 具 体 来 
说 ,中 国 东南 部 、 东 北部 以 及 中 部 靠 北 地 区 陆地 水 
储量 变化 在 滞后 1 个 月 时 与 SPEI 相 关 性 最 高 ,云南 
省 西北 部 以 及 青藏 高 原 东 南部 陆地 水 储量 变化 在 
滞后 3 个 月 时 与 SPEI 相 关 性 最 高 。 

由 表 6 可 知 ,在 中 国 陆地 水 储量 变化 与 NDVI 存 
在 0 个 月 滞后 的 顶 格 数目 最 多 , 占 比 达 到 42.0% ; 与 
降水 量 存在 1 个 月 滞后 的 栅 格 数目 最 多 , 占 比 达到 
42.1%; 与 温度 存在 0 个 月 滞后 的 栅 格 数 日 最 多 , 占 
比 达到 43.7% ; 与 SPEI 存 在 1 个 月 滞后 的 栅 格 数目 
最 多 , 占 比 达到 30.8%。 


表 6 不 同 滞后 时 长 下 的 栅 格 数目 占 比 


Tab.6 Proportion of grid number under different 


lag length /% 
影响 区 不 同 滞后 时 长 下 的 栅 格 数目 占 比 
0 个 月 滞后 “1 个 月 滞后 2 个 月 滞后 “3 个 月 滞后 
NDVI 42.0 40.0 9.2 8.8 
降水 量 24.3 42.1 23.0 10.6 
温度 43.7 9.0 18.2 29.2 


SPEI 25.3 30.8 15.4 28.5 


2.4 3 种 重要 性 分 析 结 果 比 较 对 比 地 理 探 测 器 、 
随机 森林 和 Pearson 相关 分 析 3 种 方法 对 各 个 影响 
因子 对 陆地 水 储量 变化 影响 程度 大 小 的 排序 ( 表 1、 
表 4 和 表 5) 可 知 ,地 理 探测 器 和 随机 森林 得 到 的 结 
果 总 体 相 差 不 大 ,但 和 Pearson 相关 分 析 结 果 相 比 相 
差 较 大 ,原因 在 于 Pearson 相关 分 析 只 能 反映 2 个 变 
量 间 线性 相关 程度 的 强 弱 ,而 随机 森林 和 地 理 探测 
圳 并 不 局 限于 线性 关系 。 结 合 3 种 方法 得 到 的 结 
可 知 ,气象 因素 对 陆地 水 储量 变化 的 影响 最 大 ,下 
垫 面 因素 对 陆地 水 储量 变化 的 影响 最 小 ,在 气象 因 
素 中 降水 量 对 陆地 水 储量 变化 最 为 重要 。 


3 讨论 


本 文 从 气象 .社会 经 济 和 下 垫 面 3 个 方面 人 手 ， 
定量 分 析 了 陆地 水 储量 变化 的 影响 因素 ,研究 发 现 


高 于 降水 量 ”, 且 新 疆 的 农业 绝 大 部 分 为 灌溉 农 
业 , 农 业 灌 溉 在 各 种 人 类 用 水 中 占 比 最 高 ,因此 陆 
地 水 储量 整体 仍然 呈 下 降 趋势 ;类 似 的 ,华北 地 区 
虽 降 水 量变 化 不 明显 , 因 工 业 、 农 业 大 量 开 采 地 下 
水 ,导致 陆地 水 储量 不 断 减 少 ”。 此 外 ,经 分 析 得 
到 黄土 高 原 地 区 陆地 水 储量 变化 呈 显 著 减 少 趋势 ， 
但 相关 研究 表明 黄土 高 原 近年 来 植被 增加 显著 , 主 
要 原因 是 退耕 还 林 还 草 工程 实施 以 来 黄土 高 原 深 
层 土壤 储 水 量 严重 亏损 ” 。 

本 文 结论 是 以 中 国 陆地 地 区 作为 整体 进行 影 
啊 因 素 分 析 所 得 到 的 ,对 于 各 个 地 区 , 因 气 候 不 同 、 
经 济 社会 发 展 程度 不 同 .地 理 因 素 不 同 ,各 个 影响 
因子 以 及 因子 的 交互 作用 对 各 地 的 影响 不 同 ,后 续 
人 研究 将 进一步 分 析 提 取 中 国 陆地 水 储量 变化 热点 
区 ,以 热点 区 为 单位 探究 同一 因素 对 不 同 地 区 的 作 
用 效果 ,深层 次 研究 各 热点 区 引起 陆地 水 储量 变化 
的 主导 因素 。 在 选择 影响 因素 时 , 因 可 获取 数据 有 
限 , 未 考虑 生态 输 水 产生 的 影响 中 以 及 工业 用 水 、 
农业 用 水 等 ”的 影响 ,如 何 科 学 补充 影响 因子 ,使 
之 更 具 完 整 性 , 仍 需 进一步 考虑 。 


4 结论 


(1) 中 国 大 部 分 地 区 陆地 水 储量 变化 明显 ,其 
中 显著 增加 的 区 域 主要 分 布 在 松花 江 嫩江 和 松 嫩 
平原 附近 以 及 柴 达 木 盆地 -长 江 -东南 沿海 条 带 上 ， 
显著 减少 的 区 域 主要 分 布 在 西南 以 及 新 疆 - 黄 土 高 
原 - 华 北平 原 条 带 上 。 

(2) EOF 分 析 结 果 表 明 ,中 国 陆地 水 储量 变化 
最 主要 特征 是 随 纬 度 带 交替 反 相 变化 ,其 次 是 南北 
方 陆地 水 储量 变化 相反 。 

(3) 气象 因素 比 社会 经 济 因素 和 下 垫 面 因素 对 
中 国 陆 地 水 储量 变化 的 影响 大 ,因子 两 两 交互 作用 
相 较 于 单个 因子 对 于 陆地 水 储量 变化 的 解释 力 均 
为 增加 。 

(4) 降水 量 .温度 .SPEI 和 NDVI 对 中 国 月 陆地 
水 储量 变化 的 影响 存在 滞后 性 , 且 在 不 同 地 区 各 因 
子 对 月 陆地 水 储量 变化 影响 的 滞后 时 间 各 不 相同 。 
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Changes and influencing factors of terrestrial water storage in 
China based on GRACE satellite data 
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Abstract: Determining the spatial distribution characteristics and changes in terrestrial water storage and under- 
standing the reasons behind these terrestrial water storage changes (TWSC) are necessary for the sustainable and 
comprehensive management of water resources. Based on the data of the TWSC obtained by the gravity recovery 
and climate experiment satellite retrieval, this study first analyzes the trend and spatiotemporal variation charac- 
teristics of the TWSC in China using the Mann-Kendall trend test and empirical orthogonal function (EOF) analy- 
sis. Subsequently, 10 influencing factors were selected to comprehensively analyze their relationship with the 
TWSC by employing the following three methods: geographic detector, Pearson correlation analysis, and random 
forest. The 10 influencing factors were temperature, precipitation, standardized precipitation evapotranspiration 
index (SPEI), area proportion of impervious layer, area proportion of water body, normalized difference vegeta- 
tion index (NDVI), elevation, slope, gross domestic product (GDP), and population. The results showed that areas 
with a significant increase in terrestrial water storage were mainly distributed in the areas near the Songhua River, 
Nenjiang River, and Songnen Plain, and the belt of the Qaidam Basin- Yangtze River-southeast coastal region, 
while areas with a significant decrease in terrestrial water storage were mainly distributed in southwest China and 
the belt of the Xinjiang-Loess Plateau-North China Plain. From high to low latitudes, the terrestrial water storage 
showed an alternating change pattern of high-low-high-low. Overall, meteorological factors had the strongest ex- 
planatory power for the TWSC, followed by socioeconomic factors and geomorphologic and geologic factors. 
Lag-correlation analyses showed that the monthly TWSC had a time lag response to precipitation, temperature, 
SPEI, and NDVI. The time lag of the monthly TWSC for each factor was different in the different regions. The re- 
sponse of TWSC to precipitation and SPEI mainly showed one-month lag, and the response of TWSC to tempera- 
ture and NDVI mainly showed no lag (1.e. 0-month lag). 


Key words: terrestrial water storage change; GRACE satellite; impact analysis; time-lag effect 


